A simple semiclassical approach to Kramers' problem 



Jyotipratim Ray Chaudhuri, Bidhan Chandra Bag and Deb Shankar Ray 
Indian Association for the Cultivation of Science, Jadavpur, Calcutta 700 032, INDIA. 

Abstract 

We show that Wigner-Leggett-Caldeira equation for Wigner phase space dis- 
tribution function which describes the quantum Brownian motion of a particle 
in a force field in a high temerature, Ohmic environment can be identified as 
a semiclassical version of Kramers' equation. Based on an expansion in pow- 
ers of % we solve this equation to calculate the semiclassical correction to 
Kramers' rate. 
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I. INTRODUCTION 



The recent two decades have witnessed a phenomenal development of the theory of noise 
induced escape from metastable states. The escape is governed by Brownian motion in 
addition to the characteristic dynamical motion of the system in presence of the potential 
field, V(x), Brownian motion being driven by thermal forces which, in turn, are associated 
with the dissipation 7 through the fluctuation-dissipation theorem at finite temperature 
T. The problem.0 has been addressed by a large number of workers both from classical 
and quantum mechanical point of view at various levels of description under a variety of 
reasonable approximations in the weak, strong and intermediate damping limits. The rate 
theory thus constitutes a well-developed body of literature by now!!. 

Although a number of interesting approaches to quantum theory of dissipative rate 
processes based on dynamical semigroup method for evolution of density operator were 
proposed!'! in seventies to treat the quantum and nonlinear optical phenomena with con- 
siderable success, the method could not gain good ground in the theory of rate processes 
due to the fact that it is based on weak system- reservoir coupling (7 < w). The method 
which got major appreciation afterwards in eighties and nineties in the wide community 
of physicists and chemists is the real time functional integrals. This method has shown to 
be most effective in treating quantum transition states^!, dissipative coherence quantum 
effects as well as the incoherent quantum tunneling processes!'!!!^. The latter processes 
being more relevant in the present context. Notwithstanding its phenomenal success, it may, 
however, be noted that compared to classical Kramers' theory the method of functional in- 
tegrals for calculation of rate rests on a fundamentally different footing. While the classical 
theory is based on the differential equation for evolution of probability distribution function 
of a particle executing Brownian motion in a force field, path integral methods rely on the 
evaluation of quantum partition function of the system interacting with a bath of harmonic 
oscillators with a characteristic distribution of frequencies^. The question is whether there 
is any natural extension of the classical method to the quantum domain or in other words 
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what is the analogue of Kramers' equation in the quantum regime. For this one has to look 
for an equation for evolution of a quantum probability distribution function of a particle in 
a force field driven by quantum Brownian motion. Unfortunately the theory of quantum 
Brownian motion in terms of phase space probability functionQ is far from complete. Under 
restricted conditions one may use either weak damping limit (7 <C u) or take resort to 
semiclassical procedure (ksT ^> Tiuo) and arbitrary damping. Thus, strictly speaking, to the 
best of our knowledge, there exist no full quantum analogue of Kramers' equation valid for 
arbitrary damping and temperature. 
Our aim in this paper is two-fold : 

i) Although there exists no full quantum analogue of Kramers' equation, we enquire 
whether there is any semiclassical Kramers' equation. We show that the Wigner-Leggett- 
Caldeira [WLC] equation [ see Eq. (|2.3j ) ] for Wigner probability phase space function which 
describes the quantum Brownian motion of the particle in a force field in a high temperature 
Ohmic environment may be interpreted as a semiclassical Kramers' equation. WLC equation 
reduces to classical Kramers' equation in the limit % — > 0. 

ii) While the existing methods of calculation of quantum Kramers' rate are based on 
path integral technique, we solve the semiclassical Kramers' equation for barrier crossing 
dynamics as a boundary value problem to calculate the leading order (fo 2 -order) quantum 
correction to classical rate. WLC equation had been used earlier in connection with quantum 
decoherence0 and other problems0'0. To the best of our knowledge the implementation of 
a differential equation-based approach for the calculation of quantum Kramers' rate like the 
present one had not been tried upto date. 

The outlay of the paper is as follows. In Sec. II we discuss the quantum dynamics of a 
dissipative system and show that a WLC equation for Wigner probability density function 
may serve as a good description for noise-induced escape rate in the semiclassical regime. In 
Sec. Ill elementary formulas are summarized. We apply the appropriate boundary conditions 
and make use of an ^-expansion to solve the differential equation for stationary probability 
function in Sec. IV. In Sec. V we calculate the steady state flux over the barrier, population 
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inside the well and the escape rate and show that the quantum contribution appears in the 
order h 2 . The paper is concluded in Sec. VI. 



II. QUANTUM DYNAMICS OF A DISSIPATIVE SYSTEM 

We consider a dynamical system characterized by a potential V(x) coupled to an en- 
vironment. The latter is modeled as a set of harmonic oscillators with frequencies {fij}. 
Evolution of such an open quantum system has been studied over the last several decades 
under a variety of reasonable assumptionsiHlEil. Specifically interesting is the semiclassi- 
cal limit of a high temperature Ohmic environment0. The dissipative time evolution of the 
Wigner distribution function W(x, v, t) for the system with unit mass (m = 1) for arbitrarily 
strong damping in this case can be described byQ 

(dW\ d(vW) B 2 W , . 

{-wj, , - 2 ^ +D ^ • < 2J > 

V / dissipative 

where 7 and D are the dissipation constant and the diffusion coefficient, respectively. The 

microscopic structure of 7 is given in Ref.0. x and v are c-number co-ordinate and velocity 

variables. The first term in Eq. fl2.lp is a direct consequence of the existence of 7-dependent 

term in the imaginary part of the exponent in the expression for the propagator for the 

density operator in the Feynman- Vernon theory and has been shown to be responsible for 

appearance of a damping force in the classical equation of motion for the Brownian particle to 

ensure quantum-classical correspondence. 7 and D are related by the fluctuation-dissipation 

relation, D = jhu coth 2 ^ T (in the semiclassical limit D — > 2-yksT and W(x,v,t) reduces 

to classical phase space distribution function), u is the renormalized linear frequency of the 

nonlinear system. The quantum dynamics of the bare system on the other hand is governed 

by the celebrated Wigner equation^ as given by 

(dW\ _ dW_ &VdW_ v h 2n (-l) n d 2n+1 V d 2n+1 W 

\dt ) dynamical ~ ~ V ~dx~ + ~dx~lh + ^ 2 2 "(2n + 1)! dx 2n+l dv 2n+1 ' ^ ' ' 

The above c-number description ( |2.2p is equivalent to full Schrodinger equation. The 

quantum correction to classical Liouville motion is contained in the ^-containing terms in 



the sum. We emphasize here that while ( |2.1| ) is a semiclassical Brownian dynamics due to 
the coupling of an Ohmic environment to the system, Eq ( |2.2| ) takes into account of the full 
quantum nature of the system. The overall dynamics is then given by 

(dW\ (dW\ (dW\ 



dt J y dt J dynamical \ / 



dynamical \ / dissipative 

or 

dW dW dV dW v ft 2n (-l) w d 2w+1 \/ d 2 " +1 W/ g(^W] 

~ + ~dx~dv + ^ 2 2 ™(2n + 1)! dx 2n +! dv 2n + l + 7 dv + ' ( ' ^ 

The WLC equation, ( |2.3|) is the starting point of our further analysis. Before proceeding 



further several points are to be noted : 

First, strictly speaking 7 and D terms in Eq.( [2.1|) are valid in ( |2.3| ) if the system operators 
pertain to a harmonic oscillator. When the system is nonlinear, as in the present study (and 
in the overwhelming majority of the cases in quantum and nonlinear optics), the usual 
practice is to assume that the dissipative terms remain unaffected by the Wigner higher 
derivative terms in Eq. ( |2.2| ). The validity of this assumption was examined earlier by Haake 
et. al.0 and also by us0. It is now well-known that the assumption is quite satisfactory in 
semiclassical limit and when nonlinearity is not too strong. Since the calculation of barrier 
crossing dynamics needs the linearization of the potential at the barrier top as an essential 
step, the nonlinearity of the potential V(x) must be weak. This assumption is essential for 
validity of Eq. (|2.3| ). We also mention here that weakness of nonlinearity is essential for the 
linearized description of the potential near the extrema in almost all theories including path 
integral techniques upto dateEH^I. 

Second, the system-reservoir dynamics as governed by Eq. (|2~3|) or ( |2.1| ) is based on 
Markov approximation which implies that the underlying stochastic process due to reservoir 
is Markovian by construction as in classical theory. However in deriving Eq.( p.3| ) or ( |2.1| ) the 
system-reservoir coupling is assumed to be arbitrary^. Thus Eq. (|2.3j ) is applicable to both 
in the strong and weak damping limit. Since we are considering here the spatial-diffusion 
limited rate we restrict ourselves from intermediate to strong damping range. 



Eq.( |2.3| ) has been used earlier in several occasions. For example, Zurek and have 
studied some interesting aspects of quantum-classical correspondence in relation to deco- 
herence. Based on this equation a theory of quantum noise in classically chaotic dissipative 
systems has been formulated0S The primary reason for choosing Eq.( |2.3|) as our starting 
point is that it reaches the correct classical limit, i.e., the classical Kramers' equation when 
h — > and D becomes a thermal diffusion coefficient in the high temperature limit and 
W(x,v,t) reduces to classical probability distribution function in phase space. Thus the 
WLC equation (|2.3|) may be interpreted as a semiclassical Kramers' equation. 

Our object here is to search for a semiclassical solution in powers of h for the WLC 
equation fl2.3|) under appropriate boundary conditions and calculation of the barrier crossing 
rate in the intermediate to strong damping regime. 



III. THE MODEL, THE STATIONARY FLUX AND THE WELL POPULATION 

We now consider a particle of unit mass moving in a cubic potential V(x) of the form 

V(x) = -~Ax 3 + Bx 2 , (3.1) 
3 

where x now corresponds to the reaction co-ordinate and its values at the extrema of the 
potential at x = and x = Xq denote the reactant and the transition state, respectively. In 
this model all the remaining relevant degrees of freedom of the system and the environment 
constitute a heat bath at a finite temperature T. Eq.(|2.3|) with Eq. (|3.1|) provides a complete 
description of the stochastic process in terms of Wigner probability density function of the 
c-number variables x, v. 

One of the decisive advantage of using Wigner's phase space language is that it offers 
an excellent opportunity to extend the classical formulae to the quantum domain when one 
replaces the classical probability functions by the Wigner density functions. Thus the origi- 
nal reasoning of Farkas and Kramer sSi can be followed in the present problem to determine 
the steady state escape rate from the well by considering a steady state situation in which 
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a steady state probability current is maintained between the source and the sink, the latter 
being situated beyond the transition state. The standard expression for escape rate in terms 
of flux over the population is thus given by 

k = — (3.2) 

where 

/+oo 
dvvW(x ,v) (3.3) 
-oo 

and 

n a = dx dvW(x,v) (3.4) 

J left well J 

and W(x,v) is the steady state Wigner probability density function. 

IV. THE STEADY STATE SEMICLASSICAL DESCRIPTION : BOUNDARY 

CONDITIONS 

The search for a semiclassical probability distribution function W(x,v) for Eq.( |2.3|) es- 
sentially rests on developing W(x,v) in a power series of h. Thus we writeEl 

W(x, v) = e~ Pt [p (x, v) + h Pl (x, v) + h 2 p 2 (x, v) + . . .] (4.1) 

where (3 = 1/ksT and e = ^v 2 + V(x). As pointed out by Wigne p k will not involve 
higher derivatives of V than k-th, nor higher powers of v than k-th. This is relevant for the 
explicit form of V(x) as given by Eq. (|3.1| ) in the present case and for calculation of higher 
terms in Eq. (|4.1|) in a simpler way. 

By inserting ( |4.1| ) in Eq.( |2.3| ) and equating the coefficients of different powers of h in the 



usual way and putting ^ = we obtain the following equations, 



D d 2 p | \dV \ dp y dpo =Q 
dv 2 I dx I dv dx 
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dv 2 dx 



2^v 



dpi _ n dpi 
dv 



dx 



(4.2b) 



D** + 

dv 2 



\ dx j dv 



dp 2 _ 1 d 3 V 
4! 'dx 1 ' 

f 127 2 



9a; 



12f 
Z> 2 



§7 3 3 

nV 

D 3 



A) 



, .._ , 2 67 1 dp 67 <9> . d 3 p 

+ < V > V 

\ D 2 D dv D dv 2 dv 3 



(4.2c) 



and so on. 

We note that the Eqs.( |4.2aD ,( |4.2b| ) are homogenous partial differential equations while 
( |4.2c| ) is an inhomogenous one. The inhomogenous contribution begins to appear from 
P2-term which is associated as a coefficient of h 2 in the expansion (|4. 1|) . The leading or- 
der quantum contribution (since p\ does not contribute@) through p 2 involves only cubic 



derivative of the potential in ( |4.2c ). The other derivatives of any arbitrary potential V(x) 
contribute to the terms containing higher powers of h (e.g., h 4 , h 6 etc.). Thus the leading 
order nonlinearity that affects the dynamics at the semiclassical level must at best be cubic 
one. 

Eq.( |4.2a| ) is the evolution of familiar dynamical contribution from classical Kramers' 
equation. The leading order quantum correction starts at fr 2 @. Thus W has the form 



W(x,v) 



p (x,v) + h 2 p 2 (x,v) + 



(4.3) 



where p and p 2 are the solutions of the equations (|4.2a|) and ( |4.2c|) respectively. 

To obtain the steady state probability density W(x,v) several requirements0 must be 
fulfilled: 

(i) Near the bottom of the well all the particles are thermalized and are in a state of 
thermal equilibrium. In the classical case the probability density is given by the Boltzmann 
distribution. In the present semiclassical case the corresponding distribution function con- 
tains an appropriate Wigner quantum corrections of h 2 order to Boltzmann distribution. 
According to Wigner the semiclassical equilibrium distribution is given by 



W BW (x,v)=e-^ + hf 1 + ?i 2 f 2 + 



(4.4) 



which is to be obtained solution of Eq. (|2.2|) . 
Following Wigneilll one can easily show that 

/i = 

and 



(4.5) 



h = 4vr 2 e 



(5 2 d 2 V fdV' 
~ 8 dx 2 + 28 V dx . 



+ 



(3 3 v 2 d 2 V 
24 dx 2 



(4.6) 



We emphasize the distinguishing feature of the two expansions ( |4.3| ) and (|4.4]) . While 
refers to an equilibrium thermal Boltzmann-Wigner distribution, ( [4.3[ ) denotes a steady 
state Wigner-Kramers distribution which takes into account of semiclassical Brownian dy- 
namics. ( |4.1|) may be viewed as a generalization of Kramers' ansatz in the quantum me- 
chanical context. 

We make use of the solutions ( |4.4| - |4.5|) near the bottom of the well. We thus impose the 
following condition on W(x, v) of ( |4.3| ) as 

W(x,v) — > Wbw(x, v) as x ~ and all v . (4.7a) 

This, in turn, implies po(x,v) — > 1 and p2(x,v) — > e + ^ e /2- 

(ii) On the other hand beyond the transition state, i.e., x = x all the particles are 
removed by the sink, or, 



W(x,v)xiO as x > Xq for all v . 



(4.7b) 



Having specified the boundary conditions ( 4.7a , 4~7E ) for Wigner distribution function 
W(x,v) we now elucidate the dynamics near the barrier at x — xq. Near this point we 
assume the usual linearized potential (with frequency u b ), 

1 



V(x) 



where E is the barrier height and 



(4.8) 



u 2 b = -V"(x ) > 



(4.9) 



The probability distribution function W(x,v) is determined by po(x,v) and p 2 (x,v) which 
obey the following stationary equations around the barrier ; 



D-^ + {-u 2 {x - x ) - 2-fv} ^ - = 



dv 2 



dpo dp 
dv dx 



(4.10) 



and 



D 



d 2 p 2 

dv 2 

-1.4 
12 



+ 



{-u b (x-x )-2 1 v\ — -v— 



127 



<9t> <9x 
12 7 2 



/)2 U ~ £^3 ^ \P0 +\ - 6 7 ) — - 67^ 



dpo 
dv 



d 2 po D d 3 p 



dv 2 



dv 3 



(4.11) 



The inhomogenous term in Eq. (fOT ) requires the solution po(x, v ), which can be obtained 
by solving ( 4.10|) directly using Kramers' methodSS and is given by 

A 



Po(u) = C / exp 
Jo 



2D 



z 2 \ dz + C' 



(4.12) 



where Co and C' n are constants to be determined. Here 



u = v + a(x — x ) 



(4.13) 



where a is a solution of a 2 + 27a — u 2 = and A = —(27 + a). The relevant root of a which 



makes A > is a_ = —7 — yl 2 + wf), so that ^ = — 7 + v7 2 + u l- The requirement that 



Po 



as x 



OG 



(4.14) 



implies Cq = Coi/^T - Thus 



Po(u) = C 



2A~ 



exp 



2; > cfe 



(4.15) 



Eq. ( |4. 15| ) is essentially the Kramers' solution which we now use to determine the inhomoge- 
nous terms in equation for p 2 (Eq.( |4.1ip ). Three pertinent points are to noted here. Since 
we are working in the semiclassical limit k^T > fuv, we keep the terms upto the O (rj~f) 
in the inhomogenous part. Second, we take into account of the Wigner's remark^! on quan- 
tum contribution contained in p 2 term. The term has been "interpreted as meaning that a 
quick variation of the probability function with co-ordinates is unlikely, as it would mean 
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a quick variation, a short wavelength, in the wave functions. This however would have the 
consequence of high kinetic energy" . The finite range of co-ordinates therefore ranges over 
r^Tijv where v is mean momentum ~ (/cbT) 1 / 2 . Thirdly, we also note that the leading order 
nonlinearity enters through the cubic derivative of the potential (|3.1| ) in the inhomogenous 
term. 

Based on these consideration we are now led to the following equation for p 2 : 



d 2 p 2 A dp 2 , , , „ , 9 . / Am 2 " 

+ —u-^— = ( -Mi - M 2 u - M 3 u 2 ) exp ' 



dv? D du 



2D 



(4.16) 



where the transformation of derivatives has been carried using the change of variables 



through the linear combination of x and v as expressed in ( 4.13[) . Mi, M 2 and M 3 are 
given by 

•< 2 - A 1 



M 2 = AC 



and M-t = A Cn 



7A 
A 2 



(k B T) 



1/2 



(4.17a) 
(4.17b) 
(4.17c) 



It may be noted that all the three terms in the source term of Eq. ( |4.16| ) are within (l/fe^T) 2 
order (This may be checked by noting that u behaves as v near x = xq and v ~ {kBT) 1 / 2 as 
argued earlier). 

The general solution of Eq.( |4.16| ) subject to boundary condition, 



f>2 



as x 



00 



(4.18) 



is given by 



P2{U) 



M\D M 2 D M 3 D 9 
H — ~r — u H — ~ — u 



A 2A 
M 2 D 



3A 



+ c 2 - 



2A 



exp 



2D 



2M 3 D 2 
3 A 2 

I dz + \l 



exp 



2A 



Xu 2 



C 2 - 



2D J 

M 2 D 



2A 



(4.19) 



where C 2 is an integration constant to be determined and u = v — |a_|(x — xo) with |a_ 
7 + -J7 2 + lu 2 . Thus upto h 2 our distribution function (^4.3|) becomes 
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W{u) = e"^ [ p {u) + h 2 p 2 {u) ] , (4.20) 

where po(u) and P2(u) are given by (|4.15| ) and ( |4.19|) , respectively. 

To determine the two unknown constants Co and C 2 we now employ the boundary 
condition ( |4.7a|) , which implies that the probability density function W(u) of (|4 . 2 0|) must 
coincide with the thermal Boltzmann-Wigner distribution Wbw{ x , v ) of Eq. (|4.4j) near x = 0. 

Thus we have from ( |4.20| ) [i.e., from (|4.15| ) and ([4.19|) 1 



X— > — OO _ Re 

— > e p 



C \ 



h 2 (c 9 



(4.21) 



W(x,v) , ~ u y A ■ - \- 2A J V A 

on the other hand Wbw( x , v ) of ( |4.4| ) [ with (|4.5|) and ( [4.6| )1 near the bottom of the well at 
.t = where the linearized potential 

V(x) = ~u, V (4.22) 

may be used, behaves as the Boltzmann-Wigner distribution of the particles at thermal 
equilibrium ; 

' ' (4.23) 



T T7 / \ Near x=0 _fl f 

W BW (x,v) — > e " e 



1. 
—i 

3 



1 — — II 71 p CUq 



Matching of flOID and (ggg ) suggests that 

Co = ] J 

and 

^2 



A 



2nD 



A 



^2/o2, ,2 



(4.24a) 



(4.24b) 



2A 7 3V 2nD 

The boundary condition ( |4.7a| ) therefore uniquely determines the probability distribution 
function (|4T20| ) where p and p 2 are given by ( P~T5p and (pTT9| ) and C and C 2 by ( ^24a| , ^2lrj|) 
respectively. The resultant unnormalized expression is as follows: 

( A 



W(x, v) 



+ti 2 



1 



' A 
2txD 



'ttD 
2A 



exp 



D 



dz 



cxp 



A 

D' 



-u 



M 1 D , M 2 L> , M 3j D 2 ,2 M 3j D 2 



A 



+ 



2A 



u + 



3A U + 3 A 2 



A 



3 V 2-nD 
with u = v — |a_|(x — Xq). 



it 2 (3 2 uol J exp 



.JlJ] dz --tc 2 (3 2 u 2 
2D 6 1 



(4.25) 
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V. THE ESCAPE RATE 



In order to determine the semiclassical rate of escape we now make use of the flux- 
over-population formula (|3.2| ). To this end we obtain the well population from Eqs.( p.4|) 
and (|4.25| ), assuming the linearized potential V(x) = \oJqX 2 near the bottom of the well at 
x = 0. 



n n 



2vr 

UJqP 



1 _ -K 2 tt 2 2 u: 2 
3 1 



(5.1) 



The probability current ( |3.3|) becomes 

-/3E 



J 



X 1 



7T 



+ h 2 



2\ ± + pue- 



71 



+ - A - 

2 ' 2D 



1 

3 



^2 , ,2 n2 1 

2~kD ~f3 



7T 



2 2D 



(5.2) 



where we have assumed the linearized potential V(x) = E — \u) 2 {x — xo) 2 at the barrier top 
in the expression for probability density function (|4.25|) . 

It is thus interesting to note that the quantum correction to classical population and flux 
at the semiclassical level (7i 2 -order). These yield the expression for rate of escape ( |3.2| ) as 

.4 



k 



-E /k B T 



27T UJf, 



l 



(5.3) 



where we have put (3 = l/k B T } X = —7 + y 7 2 + u 2 and A is the nonlinear parameter of the 
cubic potential ( [3.1|) . 

The above expression can be rewritten as 



k — k c i + k s 



(5.4) 



where k c \ and k sem i correspond to classical and semiclassical Kramers' rate, respectively, 
with 



^0 + e - Eo/kBT 



2tt 

UJQ 



Ah 2 



7 2 + ~ 7 



-E /k B T 



27TUb V 8 (fc s T) 3 / 2 



7 2 + ^ 2 + 7, 



(5.5a) 
(5.5b) 
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The result ( fT3| ) describes the spatial diffusion controlled rate of escape at moderate to strong 
friction 7 at the semiclassical level where the leading order quantum correction makes its 
presence felt. We also emphasize that since we are using high temperature (ksT > tkv) 
semiclassical surrounding due to Caldeira and Leggett, the result cannot be extrapolated to 
zero temperature/low temperature limit, when non-Markovian effects also begin to influence 
the dynamics. The temperature dependence of ^-containing term in (|5.3|) signifies that as 
ksT becomes large the quantum contribution significantly diminishes - a typical semiclassical 
feature in the dynamics, as expected. For 7 = the expression (|5.3|) yields the result 
of "simplest" transition state theory &tst which includes the quantum correction at the 
semiclassical as well 

= £ { 1 + d&w } — • 

For the strong friction case, Eq. ( p.3| ) can be simplified further to give the rate of escape 
in the over damped regime (7 ^> u>b) , 

0^6 f 1 h 2 A 1 Eo/kBT 

^ over damped ^ j ^ + 32 7 2 (fc s T) 3 / 2 J ' 1 ' 

which approaches to zero as 7 — > 00. h 2 term in the above expression reveals a stronger 
friction (I/7 2 ) dependence, which implies that quantum correction term is more susceptible 
to damping. The classical limit of (|5.7| ) is the Smoluchowski limit. Therefore the result 
may be regarded as the semiclassical Smoluchowski-Kramers limit. 



VI. CONCLUSIONS 

The main purpose of this paper is to explore whether a differential equation (for proba- 
bility distribution function) based approach is possible for calculation of quantum Kramers' 
rate. We show that Wigner-Leggett-Caldeira [WLC] equation which describes the quantum 
Brownian motion in a force field, in a high temperature Ohmic environment is suitable for 
this purpose. Based on an expansion in powers of fi we solve this equation for Wigner prob- 
ability density function of c-number variables under appropriate boundary conditions. We 
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make use of Wigner's quantum correction to Boltzmann factor to describe the thermalization 
of the particles deep inside the well. The effective dynamics is spatial-diffusion limited and 
the drift term (7) is Markovian. The validity of the theory is ensured at high temperature 
and for the intermediate to strong damping regime, where the quantum corrections upto the 
order h 2 significantly modifies the classical Kramers' results. 
We now summarize the main conclusions of the study: 

(i) We have identified WLC equation as a semiclassical Kramers' equation and calculated 
explicitly the barrier crossing dynamics. The method is closer to the classical one and is 
different from the existing path integral approaches for calculation of quantum Kramers' 
rate. While the existing procedure to calculate the total rate is to add k quant to the clas- 
sical Kramers' rate k c i, the present method yields the total semiclassical rate. In the limit 
h/ksT — > the classical rate is recovered. 

(ii) The quantum correction to classical Kramers' rate includes I/7 2 dependence in the 
overdamped limit (7 — > 00). 

(iii) The semiclassical theory can be reduced to 'simplest' quantum transition state de- 
scription for 7 = situation. The rate becomes susceptible to the ratio of nonlinear coeffi- 
cient to linear barrier frequency. 

(iv) The role of nonlinearity gets entangled with quantum correction in this description. 
This is obvious from the appearance of the higher derivatives of the potential in the Wigner 
equation. The quantum noise term appears as a source term in classical equation. 

(v) Kramers' celebrated ansatz in classical theory (the factorization of an equilibrium 
factor from a dynamical factor) is generalized in the present semiclassical context. 

(vi) Compared to other existing methods the quantum-classical correspondence in the 
problem of barrier crossing dynamics is relatively more transparent in the present approach. 

In view of the considerable development of the theory in the weak coupling limit, partic- 
ularly in the quantum optical context, we believe that the Wigner function approach like the 
present one can be suitably extended to analyze the energy diffusion rate in the semiclassical 
domain in the Markovian limit. Since at low temperatures, the non-Markovian effects begin 
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to be appreciable, one feels the need for appropriate non-Markovian quantum Brownian 
equations. In absence of these descriptions any differential equation based approach is to 
remain far from promising. 
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